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We derive the universal high momentum factorization of the fermion self-energy in the BCS-BEC 
crossover of ultracold atoms using nonperturbative quantum field theoretical methods. This prop- 
erty is then employed to compute the Tan contact as a function of interaction strength, temperature, 
chemical potential and Fermi momentum. We clarify the mechanism of the factorization from an 
analysis of the self-consistent Schwinger-Dyson equation for the fermion propagator, and compute 
the perturbative contact on the BCS and BEC sides within this framework. A Functional Renor- 
malization Group approach is then put forward, which allows to determine the contact over the 
whole crossover and, in particular, for the Unitary Fermi gas. We present numerical results from an 
implementation of the Renormalization Group equations within a basic truncation scheme. 

PACS numbers: 03.75.Ss, 05.10.Cc 



I. INTRODUCTION 

The physics of strongly correlated many-body systems 
plays a key role for the understanding of nature in a wide 
range of energy scales, occurring in situations such as 
high temperature superconductors, heavy ion collisions, 
or neutron stars. Therefore, it is desirable to identify 
generic features which arise beyond a given system un- 
der consideration, and shed light on the behavior of ap- 
parently different physical settings. The recent experi- 
mental progress in ultracold atoms research opens up the 
possibility to study quantum many-body systems under 
idealized conditions with full control over the system pa- 
rameters. 

A particularly promising observable in this context is 
the equation of state, which encodes the thermodynamic 
properties and the macroscopic phases of the many- 
body system [IH3] ■ Thermodynamic functions can be ex- 
pressed in terms of momentum space integrals involving 
single particle distribution functions. Due to the phase 
space factor V4nq 2 dq/(2irh) 3 being proportional to q 2 , 
these observables are also influenced by high energetic 
particles, and it is a relevant question to ask which quan- 
titative effect the large momentum tails may have on the 
equation of state. 

The single particle momentum distribution of a 
nonintcracting Fermi gas in equilibrium is described by 
the Fermi-Dirac distribution, which decays exponentially 
for momenta sufficiently larger than the chemical poten- 
tial. In contrast, for ultracold fermions in the BCS-BEC 
crossover one finds an algebraic decay j4j-[6] according to 

C 

= ^4 for large q. (1) 

Here n^^ a is the occupation number of fermions of 
species a and the precise meaning of "large q" will be 



specified later. It was realized by Tan [7HS] that the fac- 
tor of proportionality C, called contact, makes its appear- 
ance in several exact relations describing the quantum 
many-body system. For instance, it allows to compute 
the change of the energy density of the system due to a 
change of the scattering length according to the adiabatic 
sweep theorem 

1 dE = _g_ (2) 
Vd(-l/o) 4ttM' v ; 

Moreover, the energy density of the homogeneous system 
is given by 

E _ C ^ f^Q_^_( £A 

V 4nMa + J^J (2tt) 3 2M V M ° 1 ' 

Interestingly, there also is a connection between the con- 
tact and the shear viscosity of the system [TOj, [IT] . The 
Tan relations [T2T - [TS] have found verifications in experi- 
ments on the BCS-BEC crossover of ultracold fermionic 
atoms [T51124) . The measurement of the contact for a 
homogeneous system has recently been obtained in |25j . 

The contact C is related in a universal manner to the 
high momentum behavior of the fermion propagator. To 
elucidate this connection we consider the euclidean prop- 
agator of species a = 1, 2 for two atomic hyperfine states, 



(4) 

where Q = (qo,q) denotes both euclidean frequency and 
momentum. We set h — = 2M = 1 with fermion mass 
M and denote q 2 = \q\ 2 . After analytic continuation 
iqo — > — u>, the zeros of P^ )C \{uj q ,q 2 ) result in a classical 
dispersion relation oj q = q 2 — /i. For this reason we call 
^Vci(Q) = i<7o + q 2 — the classical inverse propagator. 
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The full fermion propagator has a nontrivial momen- 
tum dependence which goes beyond the quadratic behav- 
ior of the classical one. This is encoded in the complex 
valued self-energy E^,(Q), which we define in Eq. Q, 
and which is a function of frequency, momentum, and 
system parameters like temperature, chemical potential 
or scattering length. It is a challenging task of quantum 
many-body physics to (approximately) compute the self- 
energy for a given system under consideration. In par- 
ticular, for strongly correlated systems like the Unitary 
Fermi gas no obvious expansion scheme is at hand and 
one has to rely on nonpcrturbative methods like Monte 
Carlo simulations [2HH22- self-consistent T- matrix ap- 
proaches [H [5j I30H53] , 1 /W-expansion [34] , Schwinger- 
Dyson equations (SDE) [35] or the Functional Renor- 
malization Group (FRG) [37P0] , 

Using SDE and the FRG we derive the universal fac- 
torization property of the fermion self-energy 



4C 



5jjL for large q, 



(5) 



where q has to be larger than any of the physical scales 
fcph set by inverse scattering length a -1 , chemical poten- 
tial fj, and temperature T. The first contribution in Eq. 
((5), which we refer to as Tan term, results in the large 
momentum decay of the momen tum distribution accord- 
ing to n^fia = C/q 4 , see Eq. (26) below. The second 



term constitutes a high momentum shift of the effective 
fermion chemical potential (Hartree shift). The impor- 
tance of the double fraction structure of the fermion prop- 
agator in order to quantitatively describe the BCS-BEC 
crossover has first been pointed out by Haussmann in 
Refs. [H [5j . Universal high momentum factorization of 
dynamic quantities has also been observed with the FRG 
in the context of finite temperature Yang-Mills theory 
|41j and with the Similarity Renormalization Group ap- 
plied to iV-body systems at zero temperature, including 
deuteron, ultracold fermions and the electron gas [4l2"ll43| . 

Although the contact C appears in the high momen- 
tum propagator of the theory, it is a many-body quantity 
dominated by interaction and many-body scales. This 
is also resembled in its close relation to thermodynamic 
quantities. Thermodynamic considerations based on the 
Tan relations restrict the most general form of the contact 
as a function of temperature [33J EI] • I n order to com- 
pute the function C (//, T, a) it is therefore mandatory to 
work within a setting which can resolve the physics on all 
scales of the theory. This is realized by an FRG frame- 
work, where fluctuations on distinct scales are integrated 
out successively, and which goes beyond the mean-field 
approximation. 

A main result of this work is the development of an 
FRG scheme to compute the physics related to the con- 
tact. In particular, we derive a flow equation for the flow- 
ing contact Cfe, which interpolates between C'k=A = in 
the ultraviolet and the physical contact Ck=o — C in the 
infrared. The method readily applies to nonzero values 
of the crossover parameter (fepa) -1 and any spatial di- 



mension d. It allows for improving quantitative precision 
by using more elaborate truncations. We come back to 
this point in our discussion in Sec. [VI] For the purpose of 
the present work we restrict to d = 3. Moreover, we only 
consider the spin-balanced case of an equal population of 
fermions of species a = 1,2. 

We normalize the momentum distribution such that 
the integral n = 2 }^n$ a yields the density of atoms n, 
and thus the contact C is an intensive quantity. Defining 
N = 2 l~fiqa instead, with the particle number N, results 
in an extensive contact C — CV, where V is the volume 
of the system. We have C/Nk F = 3iT 2 C/k 4 . 

An interesting question is related to the range of ap- 
plicability of the asymptotic formula ([5| for the fermion 
propagator, i.e. the momentum scale q c such that the 
universal scaling form is valid for q > q c . This has direct 
consequences for the density of the system. Indeed, as 
was already noted in [5], the contribution from high en- 
ergetic atoms to the total atom density for T = can be 
approximated by 



Sn^ = 2 



d 3 g C 



C 
n 2 q c 



(6) 



The apparent divergence of this expression for q c — > 
is cured at any nonzero temperature T > and 5rS c ^ 
remains finite. Thus, there is no a priori lower bound for 
ge- 
lt can be understood easily that q c is smallest for the 
Unitary Fermi gas. Indeed, the universal form of the 
self-energy in Eq. ([5| is valid for q being larger than 
the physical scales fc p h- Typically, the inverse scattering 
length a -1 is much larger than the momentum scales set 
by chemical potential and temperature. In particular, 
this is valid in the perturbative regimes on the BEC and 
BCS sides, where \a\ —> 0. However, for the Unitary 
Fermi gas we have a -1 = and the first physical scale 
is set by either fi 1 / 2 or T 1 / 2 . For this reason there is a 
huge enhancement of the universal regime where Eq. (|5j) 
is valid in the unitary limit. This is seen in Sec. |IV| 
From Eq. ^ we then conclude that the contribution of 
high energetic particles on the thermodynamic functions 
is large at resonance. 

This paper is organized as follows. In Sec. |n]we intro- 
duce the microscopic model which is employed through- 
out the paper and deduce the factorization property of 
the fermion self-energy from the SDE for the fermion 
propagator in Sec. |III| We derive the expected results 
Cbec = 47rn/a and Cbcs = 4-7r 2 a 2 n 2 for the contact 
in the perturbative regimes, with the detailed calcula- 
tion presented in App 



IV we formulate 
After gen- 



|A] In Sec. 
the FRG approach to the contact physics, 
eral remarks on the FRG we derive the factorization of 
the flow of the fermion self-energy at large external mo- 
menta. The result is extended to the ordered regime with 
a nonvanishing field expectation value and we extract a 
flow equation for the flowing contact. We show that the 
contact becomes sizeable at the many-body scales of the 
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system. In Sec. [V]we estimate the range of applicability 
of the universal scaling behavior for the Unitary Fermi 
gas, compare the zero temperature result to perturbation 
theory and determine the temperature dependence of the 
contact at resonance. We draw our conclusions and dis- 
cuss the influence of the contact on the particle density 
and the critical temperature T c /Tp in Sec. VI In App. 
iBl we relate our discussion of the two-channel model to 
the contact in a purely fermionic language. In Apps. [C] 
and [D] we present the basic truncation which is employed 
here and derive the flow equation for the flowing contact. 
We recapitulate the derivation of the contact on the BEC 
and BCS sides from the perturbative equation of state at 
zero temperature in App. IE] 



II. MICROSCOPIC MODEL 

We consider a system of ultracold two-component 
fermions described by the microscopic many-body Hamil- 
tonian 



H = I d 3 x 



^(-V^ + X^Ut^i). (7) 



Here ipa\x) is the annihilation (creation) operator of a 
fermion atom in hypcrfinc state a at position x. Scatter- 
ing between atoms is modeled by a local contact interac- 
tion with coupling constant A^a < 0, which eventually 
translates into an s-wave scattering length a. This effec- 
tive description provides a solid starting point to com- 
pute observables for experiments with dilute ultracold 
alkali atoms such as 6 Li or 40 K, since therein small dis- 
tance details of the interatomic potential need not be 
resolved. We make this explicit by defining the theory 
with reference to a highest momentum scale A. Physics 
on distances smaller than A -1 cannot be deduced within 
the pointlike approximation and requires more informa- 
tion on the atomic interactions, like the effective range 
r c or the van der Waals length £ v dw- We choose A suffi- 
ciently large but ensure A -1 3> r c ,^vdw- In the spirit of 
this consideration we refine our definition of the contact 
according to 



C 



lim q nsa- 

9«A 



(8) 



Our analysis builds on a field theoretical formulation 
using functional integral techniques. For this purpose we 
define the microscopic model of Eq. ^ in terms of the 
euclidean action 

s = J x ( ]T rAdr - v 2 - M )Vv + A^A^lfeVi) ■ 

(9) 

The Grassmann valued field ij) a ( t j %) depends on the eu- 
clidean time variable r, which is restricted to a torus of 
circumference /? = 1/T. Adjusting T and the chemical 



potential /i allows for a description of the theory both at 
zero and nonzero temperature and density in a unified 
framework. We employ the abbreviation 



A" 



At 



(10) 



We work with a partially bosonized version of the ac- 
tion Q, where the contact interaction of fermion atoms 
is modeled by exchange of a bosonic di-atom state. The 
bosonic field enters the microscopic action by virtue of a 
Hubbard-Stratonovich transformation, which yields 



(11) 

with Feshbach coupling h v between fermions and bosons, 
and detuning v = 87ra. We introduce a counterterm 8v A 
for the constant part of the boson propagator. We ob- 
serve that the four- fermion interaction has dropped out, 
which requires A^a = ~h 2 a/^a with u\ = v + 8v A . In- 
deed, performing the Gaussian functional integral over 
the bosons one recovers a fermionic functional integral 
with action (pb. Whereas the four- fermion coupling gets 
strongly affected by fluctuations, eventually resulting in 
a divergence at the critical temperature, the Feshbach 
coupling can be approximated by its classical, momen- 
tum independent value. The scale dependence of A^ re- 
sides dominantly in the scale dependence of the term v 
quadratic in the bosons. 

The boson field (p(r, x) is not a dynamical degree of 
freedom in this microscopic action which is thus equiv- 
alent to Eq. (JoJ) . However, a dynamical part of the bo- 
son propagator will be generated due to quantum and 
thermal fluctuations. On the BCS side of the crossover 
tp describes Cooper pairs of atoms, whereas it refers to 
bosonic molecules on the BEC side. Bose condensation of 
pairs, i.e. a nonvanishing expectation value of the fluctu- 
ating field <p, leads to superfluidity of the system. We will 
find below a connection between the contact C and the 
number of atoms bound in (not necessarily condensed) 
bosonic pairs, in line with Tan's interpretation based on 
the many-body Schroedinger equation [7H5]. 



III. SCHWINGER— DYSON EQUATIONS 

The physics of the contact is based on a separation of 
scales. This scale dependence is most efficiently resolved 
within a Renormalization Group analysis. Nevertheless 
it is instructive to first consider a self-consistent gap 
equation or Schwinger-Dyson equation for the fermion 
propagator. From this exact equation, the factorization 
property of the self-energy at large external momenta 
is deduced easily. Moreover, we use this formalism to 



show how the perturbative results Cbec = 47m/a and 
Cbcs = 4:ir 2 a 2 n 2 arise naturally on the BEC and BCS 
sides of the crossover after a proper ultraviolet renormal- 
ization scheme has been applied. 

From Eq. (Ill we infer the inverse classical fermion 



P 
>- 



and boson propagators, respectively, to be given by 

P^Mio-.q 2 ) = i<7o + q 2 - Mi 



(12) 



We introduce an infinitesimal contribution iego with e — ¥ 
+ in the inverse boson propagator to regularize momen- 
tum integrals involving P v . With this single modifica- 
tion we will make manifest that only one single coupling, 
the detuning u, requires an ultraviolet renormalization 
in the BCS-BEC crossover, see Eqs. flA2"| and |A3| , Al- 
ternatively, one can add a counterterm oc /aVCVv in the 
fermionic part of the microscopic action and then adjust 
/a appropriately. At nonzero temperature, the variable 
qo is discretized according to the Matsubara formalism 
to qo = 2tt(u + 1/2)T (27rnT) for fermions (bosons) with 
n E Z. We write 



Q-P 

FIG. 1: Schwinger-Dyson equation for the inverse fermion 
propagator. A single line corresponds to a classical propa- 
gator and a double line denotes a full propagator. Fermions 
and bosons are represented by solid and dashed lines, respec- 
tively. In the loop integral we have one fully dressed Feshbach 
coupling and a microscopic one. The latter is momentum in- 
dependent. 



integration according to 



p2<kl P^P - Q)P V (P) 

+hlT VL> k i p ^ p -Q) p ^ p r (15) 



d 3 q 
(27p 



where the transition momentum fc tr is sufficiently larger 
than /c p h, such that for p 2 > k 2 T the p 2 -dependence of 
P V (P) can be neglected. We can then replace 



(13) 



For T = the sum is replaced by the corresponding Rie- 
mann integral. 

Due to quantum and thermal fluctuations, the clas- 
sical propagators P^, c i(Q) and P V}C i(Q) get dressed to 
yield the macroscopic propagators P^(Q) and P V (Q). 
This correction is encoded in the self-energies / V (Q) = = h 2 T / 
Pib/viQ) - Pib/ie>.c\(Q)- The Schwinger-Dyson equation v „ J v 



P V {P) -> PAPo,P 2 = kl) =: P^txipo) (16) 

in the second integral. By adding and subtracting a con- 
venient piece we can now write T,^(Q) in terms of two 



contributions, Ev, = £ 



(1) 



,(2) 



1 



with 



1 



1 



for the full inverse fermion propagator reads 



4 2) (Q) 



p2<kl P^P - Q) \P V (P) P v ,tM 
1 



v 



PAQ) = p^AQ) 



P P^P ~ Q)P v ,tr(P0) 



(17) 



P^P - Q)P V {P) 



(14) 



Both pieces are manifestly ultraviolet finite and do not 
depend on k tr for sufficiently large k 2 T ^> k 2 



according to the graph shown in Fig. [T] The corre- 
sponding momentum integral involves the fully dressed 
propagators. It is in this sense that the equation is 
self-consistent and cannot be solved in a straightforward 
manner. In Eq. (14) we approximate the full Feshbach 



We emphasize that the splitting of the self-energy in 
Eq. r 
gent 



15) enables us to sh ow that the superficially diver- 
oop- integral in Eq. ( 14 ) is indeed finite for large q 2 . 



coupling by the microscopic one, which is momentum in- 
dependent. This approximation is justified due to the 
weak effect of fluctuations on the Feshbach coupling. 

The BCS-BEC crossover across a broad Feshbach res- 
onance can be parametrized by the physical scales tem- 
perature, chemical potential and scattering length. We 
denote the highest physical momentum scale by fc p h. In 
the perturbative regimes, this corresponds to the inverse 
scattering length, whereas this is not valid close to res- 
onance where a -1 =0. Only momenta p 2 < k 2 h in the 
loop-integral in Eq. ( 14 ) can resolve the details of the 
choice of parameters. For this reason we split up the 



This nontrivial statement results from a scale argument 
only and thus does not imply any restrictions on coupling 
strength, density, or temperature. Moreover, the partic- 
ular choice of the artificially introduced momentum k tT 
does not play a role for the final result. 

We are interested in the behavior of Y>^(Q) for large 



momenta q' 
given by 



In Y^'(Q) we can then replace P^(P-Q) by -P,/, : ci(-Q), 
because the integral is restricted to momenta which are 
small in comparison to Q and we can neglected self- 
energy corrections to leading order at high momenta. 



^> k 2 r . The scale hierarchy in this case is 



k ph <C k t 



< q 2 < A 2 . 



(18) 
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This results in the factorization 



W(Q) = 



We define the contact according to 



for large q 



h 2 

^=7 r E 



i 



p 2 <kl \P<p(P) P<p,tr(po) 



(19) 



(20) 



since the numerator in Eq. ( 19 ) is seen below to result 



in the prefactor C of the 1 jq -tail of the momentum dis- 
tribution. Again, formulas ( 19 ) and ( 20 ) only rely on 

lich gives us a 
k 2 T . The fac- 



the splitting of the loop-integral at k t 



precise notion of "large q , namely q ^> 
torization is thus a generic feature of theories which can 
be described by similar Feshbach- or Yukawa-type gap 
equations. 

We can define an effective boson occupation number 

as 



1 



1 



with total boson number density 
yields 



C 



(21) 



This 



(22) 



and hence shows the close relation between the contact 
and the number of atoms bound in bosonic pairs. Note 
that the definition of n v has to be renormalization group 
invariant and thus involves a wave function renormaliza- 
tion constant, which we set to unity here for simplicity. 
The second term in E^(Q) results in a shift of the effec- 
tive chemical potential 



4 2) (Q) 



-Sfj, for large q . 



(23) 



The correction to the momentum distribution of parti- 
cles at high momenta which results from the self-energy 
in Eqs. (19 1 and (231 is found from the generally valid 
formula 



* = - r E 



i 



P^(Q) 



(24) 



Treating the self-energy perturbatively for large external 
momentum q 2 , we find 



^"tp.qa 



T E 

n 

r E 



i 



Pifi,d(Q) - Sfx 
AC 



PiJQ)P*A-QY 



(25) 



The first two terms yield Np{q 2 — fj, — S/i) with Fermi 
function Np(z) = ie z l T + This contribution van- 

ishes for large q 2 . Evaluating the Matsubara summation 



we find 



AC 



PlJQ)P^A-Q) 



= 4C 



l-N F (q 2 -fi) , N'piq 2 -^)^ 



<f>H 



,T C 



4(q2_ M )2 
for large q 2 , 



2(q 2 - fx) 



(26) 



which justifies the identification of AC in the numerator 
of the asymptotic self-energy. Higher order contributions 
to the fermion self-energy do not enter the 1 /g 4 -tail of 
the momentum distribution. 

The formulas derived in this section become particu- 
larly simple in the perturbative BEC and BCS regimes, 
because the integrals in Eq. (17) can be performed an- 
alytically. We present the calculation in App. [3] and 
restrict here on displaying the results. The large mo- 
mentum shift of the effective chemical potential vanishes 
on the BEC side of the crossover. On the BCS side, it is 
given by 



Sfi = -\^n^ t , 



(27) 



with four-fermion coupling = 8na. Since all atoms 
are bound to dimers in the BEC limit, the boson density 
in Eq. ( 22 ) equals half the particle density and we arrive 
at 



Cbec 



Attti 
a 



(28) 



The correct prefactor results from the wave function 



renormalization of the boson propagator, see Eq. ( A4 ) 



In the BCS limit, the first nonvanishing contribution to 
the contact arises at second order in perturbation theory 
in the coupling a. Inserting the Schwinger-Dyson equa- 
tion for the boson propagator (|A2[) into Eq. (22), we 



obtain a double integral over two fermion propagators, 
each resulting in a fermion density n^, j(T . We find 



Cbcs = An 2 a 2 n 2 . 



(29) 



These perturbative results derived from Schwinger- 
Dyson equations agree with the expressions found from 
the zero temperature equation of state. The correspond- 
ing calculation is recapitulated in App. IEJ 



IV. 



FUNCTIONAL RENORMALIZATION 
GROUP 



The factorization ( 19 1 relies on the scale hierarchy ( 18 ). 
The latter is also at the root of the derivation of the 
contact C in Eq. (201 by means of the splitting of the 
integration defined in Eq. (15). A method of choice 



which naturally incorporates such a momentum splitting 
and the factorization in the fermionic self-energy is the 
Functional Renormalization Group (FRG) approach to 



G 



the BCS-BEC crossover. This has been already shown in 
the FRG approach to thermal propagators in Yang-Mills 
theory [41 . In the present context it allows us to extract 
the many-body quantity C(a, /i, T) in the whole crossover 
without the need to introduce the momentum splitting 
with fc t r- It also goes beyond the universal scaling be- 
havior proposed in Eq. ([5]). In fact, the only assump- 
tion underlying the following analysis is the perturbative 
treatment of the fermion propagator at high momenta, 
whereas the bosonic degrees of freedom are treated in 
a nonperturbative fashion. The latter is in particular 
mandatory to resolve the second order nature of the su- 
perfluid phase transition. The truncation employed in 
this section is introduced in Apps. [C|and[E>} 



A. Flow equation 

We start from the microscopic action of the two- 
channel model given in Eq. (11) at a high momentum 



scale A and successively include quantum and thermal 
fluctuations on momentum scales larger than a momen- 
tum scale k. By lowering k we resolve the macroscopic 
properties of the system in a coarse graining procedure 
and, eventually, arrive at the full effective action T of 
the theory for k = 0. The effective action generates all 
correlation functions of the system. (More precisely, it is 
the generating functional of the one-particle irreducible 
Green functions, which are obtained by simple functional 
differentiation.) In between, the effective average action 
Tfe can be interpreted as the effective action for a probe 
with size k~ x . Although for k > is not a physi- 
cal observable, as it depends on the scheme in which we 
lower fc, this rough picture provides an intuitive and of- 
ten surprisingly accurate guiding principle. In summary, 
the effective average action satisfies 



S, 

r. 



(30) 
(31) 



The evolution of as k is lowered is governed by the 
exact flow equation [45 



d fc I\#, i/>] 



i STr f L_ 

2 Vlf^] 



-dkRk 



(32) 



where STr is the "supertrace" , which takes into account 
the mixed bosonic and fermionic structure of the theory, 
and Rk is the regulator related to the mentioned way how 
k is lowered. In a momentum representation, the kernel 
Rk(Q)$(Q + Q') of the operator R k needs to satisfy 



lim R k (Q) = k 2 , 

q 2 /k 2 — >0 

lim R k (Q) = 0. 



(33) 
(34) 



The same is true in the corresponding limits of |qo|/& 2 - 
If we were able to solve the functional differential equa- 



many-body Hamiltonian or computing the complete path 
integral, we would have solved the theory. In most practi- 
cal cases, this is far beyond reach and otherwise exactly 
solvable models are typically diagonalized in a simpler 
fashion than employing the renormalization group. The 
actual strength of the FRG is that it can deal with com- 
plex systems, both perturbatively and nonperturbatively, 
since the effective average action Y k provides an intu- 
itive way of organizing the inclusion of all fluctuations 
in the path integral and is thus particularly well suited 
for physically motivated approximations: We restrict the 
space of possible functionals Y k to a limited set of corre- 



lation functions, project the flow equation (32 1 onto this 



parameter set, and solve the remaining set of equations 
numerically or maybe even analytically. Due to such a 
truncation we practically never arrive at the full effective 
action T which generates all correlation functions, but 
obtain a quantity which is usually predictive enough to 
describe, for instance, the low-energy physics and ther- 
modynamics of the system. 

The application of the FRG to a broad variety of top- 
ics, ranging from high energy physics to condensed mat- 
ter, are reviewed e.g. in [4"5H5"3"] . The FRG approach to 
many-body physics with ultracold atoms is reviewed in 
l54l and the lecture notes l55l. 



B. Fermion self-energy and contact term in the 
symmetric regime 

For conceptual clarity we first isolate the contact term 
and the shift of the chemical potential from the flow equa- 
tion of the self-energy T,^{Q) — P^,(Q) — P^,a(Q) in the 
symmetric or disordered regime, where the field expec- 
tation value of the boson field is zero. The procedure 
is easily extended to the broken regime of the flow with 
a nonvanishing expectation value and the result is given 
below. Note that in general the self-energy is a 4 x 4- 
matrix. For equal population of the two hyperfme states 
a = 1,2, the most general form of the self-energy can 
be parametrized by two complex functions S^(P) and 
£i/>,an(-P), where the second one is called the anomalous 
contribution. For simplicity we neglect the anomalous 
self-energy in the following. We show in App. [C]that it 
does not contribute to the contact. 

The flow of the inverse fermion propagator, which is 
identical to the flow of the self-energy, is given in the 
symmetric phase by 



dk^"ip,k(P) 



dkR<p,k(Q) 



1 



'P^kiQ) + R<t>AQ)) 2 P*,h{Q - p ) + r tpAQ - p) 

dkRip,k{Q) 1 

(Pi,AQ) + Ri>,k{Q)) 2 p<pAQ + p) + R<pAQ + p) 

(35) 



tion (32 1, which is equivalent to diagonalizing the full We emphasize that the external momentum P is a free 
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parameter and for each P we have an individual flow 
equation. We neglect a possible scale dependence of 
the Feshbach coupling h^. The expressions P^AQ) = 
Pij>,h{Q) + "£>ip,k(Q) and PqAQ) are the full inverse prop- 
agators at scale k. We have 



p^MQ) = i^o + q 2 



//,. 



(36) 



Since the boson propagator is gapped for k ^> fc p h, 
where fc p h is a physical scale given by either the inverse 
scattering length, temperature or chemical potential, we 
effectively only have a nonvanishing contribution to Eq. 
(351 for k < fcph- Choosing a large ex ternal frequency 
> k 2 , we can use property ( 34 ) of the regulator 



p*% k 2 ph 



functions Rk(Q) to approximate 

P k {Q ±P) + R k (Q ± P) ~ P k (±P). 



(37) 



Thus, for large external momentum P the flow of the 
self-energy simplifies according to 



d k ^, k (P) 



d k Rt, k {Q) 



p^A-P) Jq (p*AQ) + R<»AQ)f 

d k R^. k {Q) 



P*,k{ p ) Jq (Pi>AQ) + R*AQ)Y 



(38) 



We now show that this results in an asymptotic self- 
energy E-i = Ew, fc = o of the form 



AC 



5fi. 



(39) 



The first term in Eq. (38) yields the contact term. In- 



deed, for large external momenta P, renormalization ef- 
fects on the fermion propagator are small and we can 
approximate P^, jk (—P) ~ P^^—P). We then find the 
P-dependent term being multiplied by an integral which 
receives contributions from the physical scales fc p h. Since 
P*I>,a(—P) is fc-independent, we can integrate the first 
term in Eq. ( 38 ) and identify the contact as being given 
by 



AC 



dkh% 

A JQ 



(I 



hi 



dkd k , 

a JqP+M + RoAQ) 



(40) 



In the second line, we introduced the formal derivative 
d k , which only acts on the /c-dependence of the regula- 
tor R k . The advantage of this rewriting is to make the 
simple one-loop structure of the equations manifest. Eq. 
(40) allows us to define a scale dependent contact C k 



according to the flow equation 



d k C k = -^-dk 



p^AQ) + R*AQ) 



(41) 



with Ca = and Ck=o = C. 

The second term in Eq. ( 38 1 contains the boson prop- 
agator evaluated for large momentum. Since the micro- 
scopic boson propagator is constant and the momentum 
dependence only builds up due to the renormalization 
group flow, this contribution to the fermion self-energy 
is independent of P and constitutes a shift of the effective 
chemical potential. The asymptotic value is then equal 
to the one evaluated for a large momentum fc tr . We con- 
clude 



6jj, = — 



dk^-dk 



1 



p*AQ) + fyAQY 



(42) 



With the effective four-fermion vertex —h^/m^ — we 
find Sp — —X^n^a-. 

C. Extension to the ordered regime 

For low enough temperatures, a nonvanishing expecta- 
tion value po,fc 01 the boson field p — (j)*(f> appears during 
the renormalization group flow. If po.fc=o — Po > 0, we 
say that the system is in its superfluid phase and some 
of the bosons have condensed. Note that, due to inter- 
actions, the condensate fraction does not coincide with 
the superfluid fraction. Above the critical temperature 
there is a region where a nonvanishing value of po,fc ap- 
pears for k > during the flow, but does not persist for 
k — > 0. We then arrive in the symmetric (normal) phase 
of the system. We may call this intermediate region the 
precondensation regime. 

Conceptually the above derivation of the flow of the 
self-energy and the asymptotic scaling with contact term 
~ 4:C/PjjiA~P) remains valid also in the presence of 
a possibly nonvanishing boson field expectation value. 
The corresponding calculation is presented in App. |D| 
We find from Eq. (D14) that the flow equation for the 
contact is given by 



d t C k 



'{VAPo.k + d t poA 



(43) 



2 d/2 



V d > 



n 



d 



S 2 '" V d + 2 



N' B (uct>A 



<p,k 



k a+ ' ( 1 - 

d-t 

N B {^4>A - ^4>,kN' B {uj^ k ) 



with 



^4>,k — 



in 



<p.k 



w 



4>,k 



2XkPo,k 



Si- 



(44) 



t = log(fc/A), v d = (2 d + 1 7r rf / 2 r(d/2))- 1 and d = 3. The 
Bose function is defined as Nb{z) = [e z l T — 1) _1 and 
N' B = dNs/dz. For a definition of the corresponding 
running couplings we refer to App. [C] Typical renormal- 
ization group flows of Ck for the unitary Fermi gas are 
shown in Fig. [2] 
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.... t = 





— T — 


T c 


- T = 


1.2T 


— T = 


2T C 



O 400 



t = log(fc/A) 




FIG. 2: RG-scale dependence of the flowing contact Ck at 
unitarity a -1 = 0. We have k = Ae* such that t = corre- 
sponds to the ultraviolet and t — > -co to the infrared. We 
observe that the contact is unaffected by fluctuations of ul- 
traviolet modes and it starts to build up on the many-body 
scales of the system, which are set here by the chemical poten- 
tial and temperature corresponding to t M = log(/i/A) = —6.9 
and tr — t^. Obviously, all curves saturate at a certain value 
of t and we can read off the physical value at k = 0. 



At zero temperature we obtain a nonvanishing value 
for the contact C. The corresponding value is found from 



(43) by setting the Bose functions to zero. We have 



d t C k \T=o = -^(ilAPo,k + d t po,k) 



V 

d + 2 



(45) 



V. RESULTS 

Within the truncation scheme put forward in App. [C] 
we can compute the contact as a function of the crossover 
parameters /i, T, and a. Moreover, the high momentum 
factorization of the self-energy can be shown explicitly by 
solving the flow equation for S.0(P) for different values 
of P. In order to translate the results expressed in terms 
of the chemical potential for the density, the equation 
of state T) has to be applied. Since the density 
(and, iteratively, the contact itself) receives substantial 
contributions from the contact term in the fermion self- 
energy, fully self-consistent results can only be obtained 
from a self-consistent treatment of the Tan term in the 
flow equations. Here we restrict to an analysis of the 
qualitative behavior of the contact and do not aim at 
quantitative precision. 

At zero temperature, the expression for the contact on 
the BEC side can be derived from the Lee-Huang- Yang 
(LHY) equation of state (E2) and the adiabatic sweep 
theorem (2 1 
where we a 



FIG. 3: Zero temperature contact on the BEC side of the 
crossover. The many-body chemical potential is defined as 
Mmb = At — £b/2 and thus a positive quantity. (See for instance 
Eq. 



(E5l in the context of a weakly interacting Bose gas.) 
The FRG treatment cap tures the Lee-Huang- Yang (LHY) 
correction, see. Eq. (E9|. Mean field theory (MF) is shown 
by the dashed curve. 



of /i m b = P — £b/2, with (negative) binding energy £b of 
the molecules. 

The result of the integration of the renormalization 
group equations at zero temperature on the BEC side 
is given in Fig. [3| We find excellent agreement with 
the prediction from LHY theory, whereas the mean field 
curve deviates substantially. From Eq. ( D37 ) it is appar- 



ent that the LHY correction, which is reproduced in the 
equation of state on the BEC side as well, is also visible 
in the contact, because both share a common flow equa- 
tion. The nontrivial renormalization of the prefactor in 
Eq. (D37) ensures the result to be beyond mean field. 



We recapitulate this derivation in App. [Ej 
so give the corresponding expression in terms 



As we approach unitarity from the BEC side, we leave 
the perturbative regime and the contact is no longer de- 
scribed by the LHY expression. The FRG result within 
the truncation of this work is shown in Fig. |4j For 
a" 1 = we find C/p 2 = 0.34 and C/k* = 0.11 at 
zero temperature. The Bertsch parameter within this 
approximation is £ = 0.55. We observe the contact to 
be a monotonous function of the crossover parameter 
(V^mbfl) -1 , or, equivalently, (fcpa) -1 . 

For negative scattering lengths we find our zero tem- 
perature results to be far below the BCS prediction. 
The reason for the failure of the present truncation is 
that the momentum dependence on the BCS side is not 
well-approximated by a derivative expansion of the bo- 
son propagator, although momentum independent ob- 
servables like the equation of state are described cor- 
rectly. This becomes transparent in the derivation of 
Cbcs = 47r 2 n 2 a 2 from the Schwinger-Dyson equation in 
Eq. (A8), where we explicitly use the momentum depen- 
dence of the bosonic self-energy Y, V (Q). In App. |B| we 
discuss how the contact on the BCS side is more easily 
derived in a purely fermionic language. 
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FIG. 4: The contact close to resonance as a function of the 
scattering length for both T = and T = T c (a, n). The la- 
bels of the axes are analogous to Fig. [3] At unitarity we 
obtain C(T = 0)//i 2 = 0.34 and C(T = T c )/fi 2 = 0.24. Far 
on the BCS side our present truncation becomes inappropri- 
ate as is discussed in App. [B] For better visibility we show 
the asymptotic BCS value up to (^/Jla)^ 1 = —0.5, which is 
already beyond the applicability of BCS theory. 



FIG. 6: The blue solid line shows the temperature depen- 
dence of the contact normalized by the chemical potential for 
the Unitary Fermi gas. We observe a decrease of C/fi 2 as we 
approach the critical temperature from below, resulting in a 
sharp dip at T c . The function is monotonic for T/T c > 1.75. 
The green, red, and orange curves correspond to different con- 
tributions to the contact and are explained in Eqs. ( 46 (— ( 48 1 
in Sec. HVl 




1.5 

T/T c 



FIG. 5: The contact of the Unitary Fermi gas normalized 
by fcp with Fermi momentum k-p = (37r 2 n) 1//3 . Relating this 
plot to an extensive contact C = CV with volume V we 
have C/Nk F = 3-K 2 C/kp, see our discussion of the normal- 
ization in the introduction. For zero temperature we obtain 
C(T = 0)/kp = 0.11 and C(T = T c )/kp = 0.10 at the criti- 
cal temperature. The Fermi momentum is not corrected due 
to the high momentum contribution to the particle number 
density. Hence, kp will in general be larger than plotted here. 



In Fig. [4] we also show the value of the critical contact 
C(T C ) in the crossover. We find the corresponding value 
always to be below the zero temperature value. The full 
temperature dependence of the contact of the unitary 
Fermi gas is shown in Figs. [5]and[6j We observe a sharp 
dip at the critical temperature, hence C(0) > C(T C ) in 
Fig. [4] Since the contact is related to a first derivative 
of the energy (or pressure) according to the adiabatic 
sweep theorem ([2]), it has to be continuous at T c as a 



result of the second order nature of the phase transition. 
We confirm this behavior in our results. The contact 
C/kp shows a maximum above T c . The necessity of a 
maximum can be derived from general thermodynamic 
considerations [44 . 

From Eq. (43) we observe that the contact receives 



contributions from different terms in the flow equation. 
These are important in distinct regimes of the system. 
To visualize this, we split up the flow of into three 



parts according to d t Ck = d t C i 



d t c k ] = -^(riAPv.k + d t po,k), 



(i) 



d t C, 



(2) 
k 



d t ci 3) with 



(46) 



o t cP = 



2 d/2 



(48) 



The only term which persists in the symmetric phase, 

(2) 

and the stages of the flow where po,fe = 0, is . There- 
fore, it is the leading contribution above T c . Both 

(3) 

and C k start to build up in the precondensation phase. 

(3) 

However, Cf, is never really large, and the contribution 
of c\}^ above T c is negligible. In the zero temperature 
limit, the term becomes most important. This can 



be understood easily from Eq. (D31I, where we identify 



this term as the contribution from condensed bosons to 
the density. 
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• ReE^ 





FIG. 7: The asymptotic approach of the contact for the uni- 
tary Fermi gas at the critical temperature. We plot the real 
part of Pip,cl(— P)'S^(P) with P^,, c i(— P) = — ipo +P° — A* and 
po — ^T. In accordance with formula |5|, the factorization 
at large momenta leads to the approach of the constant value 
4C. For very low momenta our perturbative treatment of the 
fermion propagator becomes quantitatively less accurate, but 
still the deviations do not exceed 20% even for p — s> 0. We re- 
strict the self-energy here to the diagram in Fig. [TTJwhere the 
external momentum P appears in the fermion line and which 
is responsible for the contact term. This effectively adds the 
constant 8[i to E^, for large momenta. 



FIG. 8: The fermion self-energy T,^ t k=o{P) at a -1 = and 
T = T c computed from the first diagram in Fig. |11| We eval- 
uate the function for po — nT and find the large momentum 
behavior to be a reasonable approximation for both the real 
and imaginary parts even at low momenta. The asymptotic 
form Q is shown by a dashed line. The self-energy vanishes 
for large momenta since we effectively added the constant 5fi 
for p — > oo by neglecting the second diagram in Fig. |11| 

the reliability of this finding, but is beyond the scope of 
the present work. 



We already addressed the question whether the scal- 
ing formula ([5| can be applied for a large part of the 
momenta or only yields an asymptotic, but practically 
irrelevant contribution. For this purpose we solve the 
flow equation for the self -ene rgy E^fc(P) on a grid of P- 
values according to Eq. ( D9 ) . Therein we restrict to the 



first integral, which corresponds to the diagram in Fig. 
[Tljwhere the external momentum appears in the fermion 
line. Only this diagram contributes to the contact. The 
universal regime of validity is expected to be large for the 
Unitary Fermi gas. In Figs. [7] and [8] we underline this 
statement at T = T c and a -1 = 0. 

Since the self-energy E^(P) is a complex valued func- 
tion of P = (po,p), we gain information about the high 
momentum behavior from plotting both the real and 
imaginary part at the lowest possible fermionic Matsub- 
ara frequency po — irT as a function of To see the 
asymptotic approach of the form P^,, c \{— P)T*^ — 4C we 
plot the real part of this particular combination in Fig. [7] 
The imaginary part of this product vanishes for p — > oo, 
showing that the contact indeed is real-valued. 

We plot ReE^,(P) and ImS^,(P) in Fig. [8] for the same 
set of parameters as before. We find that the scaling form 
is a good description for all momenta at the critical tem- 
perature. Although this does not come unexpected for 
a scale invariant, critical system, this behavior could be 
a relict of our perturbative treatment of the propagator 
for small momenta, where this is not necessarily a valid 
assumption. Further improvement of the truncation and 
iterative solution of the flow equation will shed light on 



VI. DISCUSSION & OUTLOOK 

In this work we derived the contact in the BCS-BEC 
crossover from a universal factorization of the large mo- 
mentum part of the fermion self-energy. The analy- 
sis is built on exact Schwinger-Dyson equations and 
an exact Functional Renormalization Group equation. 
Both methods allow for a computation of the contact 
C(/i,T, a). In the limiting perturbative cases the ex- 
pected relations Cbec = 47rn/a and Cbcs == 47r 2 n 2 a 2 
can be obtained analytically. 

The factorization ^ of the self-energy into a universal 
part and the momentum independent many-body contact 
can be shown explicitly by solving the renormalization 
group flow equation for E^, ^(P) for particular values of 
P. The regime of universal scaling turns out to be large 
for the Unitary Fermi gas. This has profound conse- 
quences for the density. Indeed, if we assume the scaling 
behavior of the self-energy to be valid for all momenta, 
we can estimate the contribution of high momentum par- 
ticles to the total density by 



J^(C, estimate) g£» 



/ d 3 p /l 
J (2tt)3V 



-2N F (j> 2 ~y) 
4(p 2 - /i) 2 

n' f (p 2 - M) 

2(p 2 - M ) 



(49) 



see Eq. (26) 



This density correction from high momentum particles 
has a direct impact on the value of the critical tempera- 
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ture. The discrepancy of the critical temperature at uni- 
tarity T c /Tp — 0.276 neglecting the density correction, 
expressed in terms of the density, is larger in compar- 
ison to other methods than the grand canonical value 
T c / fi = 0.44. Part of this mismatch is likely to be found 
in the normalization with an incomplete density where 
substantial contributions are not included. Assuming for- 
mula ( 49 1 to be valid, we find a correction of the density 
j n (C,efenW)/( n (0) + £ n (C,e 8 timate)-) _ q.29, resulting in a 

corrected value T c /Tp = 0.22. Here r?/ ' is the density 
which is obtained from the standard truncation without 
contact as it is described in App. [C] See Eq. (D31) for 
the corresponding flow equation dtnk = dtn k °\ 

The suppression of the ratio T c /Tp due to the high 
momentum tail of the particle density distribution also 
sheds light on the success in computing the critical tem- 
perature with methods which take into account the full 
momentum dependence of the propagators, e.g. Monte 
Carlo simulations or self-consistent gap equations. How- 
ever, having singled out this important contribution to 
the fermion propagator, an effective model which takes 
into account the Tan term might give reasonable results 
with less numerical effort and, furthermore, help improv- 
ing our understanding of the physics of the BCS-BEC 
crossover. 

To further improve our truncation we can feed back the 
effect of the contact on the remaining running couplings. 
Indeed, so far we only evaluated d t Ck on the solution of 
the standard truncation. However, as a correction to the 
fermion propagator, the self-energy results in a modified 
flow of, e.g., the effective potential. A different running 
of the boson mass term k results in a different phase 
diagram and hence a corrected value for T c / fj,. 

To close the present truncation of the FRG equations 
up to order P ( ~ 1 " 1 (Q)-effects in the self-energies, we have 
to also include the bosonic contact, which arises in the 
high momentum limit of the bosonic self-energy. The 
relevance of the contact corrections is most pronounced 
for the pressure of the system, which is minus the effec- 
tive potential and thus only includes a single propagator 
in the corresponding flow equation. However, we have 
seen on the example of the density that also corrections 
to first derivatives of the effective potential, i.e. den- 
sity rife, boson mass m| k and condensate density po,k, 
can be substantial. Hence a closed set-up for improve- 
ment of our truncation needs to incorporate the fermionic 
and bosonic contact as a correction to the flow of these 
couplings. This can be achieved in a perturbative and 
iterative fashion. 

Acknowledgements 

We thank T. Enss, D. Goncalves-Netto, M. Scherer, 
D. Schnoerr, and W. Zwerger for discussions and help- 
ful comments. S. D. acknowledges support by the Aus- 
trian Science Fund (FWF) through SFB FOQUS and the 
START grant Y 581-N16. I. B. acknowledges funding 
from the Graduate Academy Heidelberg. This work is 
supported by the Helmholtz Alliance HA216/EMMI. 




FIG. 9: Schwinger-Dyson equation for the inverse boson 
propagator. The notation is chosen as in Fig. [I] 



Appendix A: Contact in the perturbative regime 



To exemplify our statements in Sec. Ill we consider 



here the perturbative BEC and BCS regimes, where the 
integrals in Eqs. ( 19 ) and ( 23 ) can be performed analyt- 
ically. For weak interactions, the largest physical scale 



is given by k ph 



To study the contact in the 



whole crossover we employ the Functional Renormaliza- 
tion Group in Sec. |IV| One of the merits of this method 
is that all expressions are automatically renormalizcd. 

On the BEC and BCS sides of the crossover, the 
fermion propagator only gets weakly dressed because ei- 
ther the fermion or the boson propagator is gapped. This 
leads to a suppression of the loop-integral originating 
from the diagram shown in Fig. [l] Hence we can al- 
ways treat the fermion self-energy perturbatively in these 
regimes. We identify the large momentum shift of the 
chemical potential as 



J P n V 



1 



(Po - Qo) + P 2 - M) (f + ie Po) 



-N F (p 2 -li)-N B {v/e) 
e(iq -p 2 + n) +v 



h 2 

<p 



(Al) 



on the BCS side. The Bose function is denoted by 
Nb(z) = [e z l T — 1) _1 . In the last expression we take 
the limit s — \ + . The go-dependence of <5/i defined in 



Eq. ( 23 ) is seen to vanish. 



We first consider the BEC limit, where fj, ~ — 1/a 2 is 
large and negative. The shift of the chemical potential 
vanishes due to the suppression of t he Fermi function 
Nf(p 2 — A 1 ) m the integral in Eq. (Al). The inverse 



boson propagator is derived from the Schwinger-Dyson 
equation [56] for the bosonic self-energy 



s v (Q) = p v {Q)-p v AQ) 



i 



(A2) 



shown diagrammatically in Fig. [9] The counterterm is 
given by 



hlA 
SiyA= 4^- 



(A3) 



In the perturbative regime we can replace the full fermion 
propagators in the loop-integral by the classical ones. 
The Matsubara summation and angular integration can 
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then be evaluated analytically. We do not need the full 
expression but only note that 



P. 



(Q)=v + 



h 2 

•p 

8^ 



2 



21 
A 



H ~ Z v I ig + 



(A4) 

for |g |, q 2 -C a -2 in the BEC regime. (We used —h^/v — 
8ira and // = — 1/a 2 ) The boson propagator resembles 
particles with classical dispersion relation u q — q 2 /2M v 



and mass M, r 



2M = 1. The wave function renormal- 



ization constant is Z v — /i 2 a/32-7r. 



Inserting the boson propagator from Eq. ( A4 1 into Eq. 



( 20 1 we arrive at 
8tt r 

a Jp 2 <k 2 r 



Cbec — 



i 



i Po + p 2 /2 



8tt 



p 2 <k\ 



(N B (p 2 /2)-N B (ki) 



ip + fcj 

87rn ViC i 
a 



(A5) 



where we used fc p h 3> T 1 / 2 and n V)C j defines the number 
density of boson with classical dispersion relation from 
Eq. (21 1. Due to the equation of state n v _ c \ = n/2 on 



the BEC side of the crossover we conclude 

Aim 



BEC 



(A6) 



as expected. The perturbative expression for the contact 
starting from the equation of state at zero temperature 
is recapitulated in App. [EJ 

In the BCS regime the bosons are resonant excitations 
and n^ cl — 0. The corresponding formula for the con- 
tact is most easily derived from inserting the Schwinger- 
Dyson Eq. ( A2 1 into formula ( 20 1 . Therein, the boson 



self-energy T, v oc h can be treated perturbatively, since 
it is small in comparison to the boson gap v due to the 
small scattering length a = —h^/Snv. We then find 



C 



BCS 



h 2 

n 

Av ^ 



E„(P)-S v (Po,^) 



P 2 <kl 



iepo) 2 
1 



P 2 <kl Jk p^A k ){ v + ' 1£ Po) 



1 



1 



Pip,d(P - K) P-4>, c \(jPo - ko, (hi 
r N F {{p-k) 2 -u) 
K J P 2 <kl P^, c \{K)(v + iek ) 



K 

Av 



K 

Av 

K 

Av 



1 



A F (fc 2 - M ) 



iefc ) 

■N F {v/e) 



(A7) 

We again applied the limit e — > + . With the BCS equa- 
tion of state n^ j(T — n/2 we arrive at 

C B cs=47rVn 2 . (A8) 



e(k 2 - n) 



This agrees with Eq. (E12) 




FIG. 10: Schwinger-Dyson equation for the inverse fermion 
propagator in a purely fermionic theory. Again, a single or 
double line denotes a classical or full propagator. The empty 
(filled) dot corresponds to a classical (full) four-fermion ver- 
tex. We can combine the last two terms such that the correc- 
tion is only given by the tadpole diagram but with a momen- 
tum dependent effective vertex. 



Appendix B: Contact in the atomic phase 

Deep in the atomic phase, i.e. for high temperatures 
or on the BCS side of the crossover, bosons are excited 
resonantly. For this purpose, an ansatz P^(Q) ~ Ptt>, c \{Q) 
fails to capture the right behavior of the contact. This is, 
however, not a failure of the Yukawa model, but rather 
shows the need for keeping higher order terms in the 
bosonic self-energy. We have seen in Eq. (A8) that the 
BCS result Cbcs = 47r 2 a 2 n 2 can be obtained from the 
full boson propagator. Our truncation for solving the 
FRG equation for the fermion propagator, as it is de- 
scribed in Sec. 



IV and App. O cannot resolve this be- 



havior, since the bosonic self-energy is only incorporated 
by the wave function renormalization coefficients A^, 
and the boson mass term mi. This is made explicit in 



Eq. (C4) 



The contact in the atomic phase can be obtained 
within a purely fermionic language. Then the momen- 
tum dependence of the related Schwinger-Dyson equa- 
tion for the propagator is solely encoded in the two-loop 
diagram shown in Fig. [TUJ because the tadpole diagram 
contains the classical vertex, which is momentum inde- 
pendent. Equivalently, we can express the correction in 
terms of a tadpole diagram with full four-fermion vertex. 
A (2PI-)resummed form of the latter is given by 



A,/ 



i-A. (n(p)-n A (p) 



(Bl) 



with U{Q) being defined in Eq. (|A7} and U A (P) cho- 
sen as a suitable subtraction at large momentum. Ne- 
glecting renormalization issues, we find for large external 
momenta p 2 to second order in Xj, that 



.(P) = I 
Jc 



Ay,(Q) 

Q PAQ-P) 

1 

P^,cl(-P) 



\2 



n(Q)-n A (Q) (B2) 



Emp loying A^, = — 87ra we recover the integral in Eq. 
flA7 ) and thus C B cs = 47r 2 a 2 ?i 2 . The FRG allows to go 
beyond the bubble resummation in Eq. (B2), see e.g. 
Ref. [57]. 
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Appendix C: Truncation of the effective average 
action 



In order to solve the flow equation (32 1, we have to em- 



ploy a truncation of the effective average action Tk [</>, ip] 
which is on the one hand tractable from a computational 
point of view, but also contains the relevant physics of 
the underlying system. For this purpose we make the 
ansatz 



(Cl) 



For a detailed discussion of this nonperturbative deriva- 
tive expansion scheme see [5SH40) . The Yukawa model 
applied to many-fermion systems has also been studied 
in Ref. [37] . The effective average potential U k (p, p) 
with p = (f>*(f> is expanded in a power series around the 
physical values /?o,fc and p according to 



Uk(p, A) 

= ™l,k.p,(p~- Po,k) + 



L 4>,k,{i\l J - W,k) -1 ^ (P ~ P0,k) 2 + ■■■ 

(™\k(P ~ Po,k) +a k (p- p)(p - p 0fk ) + 



(P-po,kr + ... + 



(C2) 



We neglect higher terms in this series and only keep track 
of the running of m 2 fc , A^fc and a&. Going beyond this 
basic set of parameters yields higher quantitative preci- 
sion. Note that we neglect a possible scale dependence 
of the Yukawa coupling Ka,. The initial conditions for the 
running couplings are set by T\ = S to be 



Z A = A A = 1, m\ A = v A - 2p, 
A^a = 0, «a = -2, p Q .\ = 0. 



(C3) 



The vacuum flow of k allows for extracting the scat- 
tering length a which is related to the detuning from res- 
onance v(B). Note that we define the detuning slightly 
different here than in Eq. (Ill since it does not contain 
the chemical potential. 

The ansatz in Eq. \G1\ takes into account the impor- 



tant low momentum modes of the boson propagator. To 
see this more clearly we expand the bosonic self-energy 
in the exact relation for the boson propagator 

P<t>{Q) = P<p,a(Q) + ^4>(lo, q 2 ) 

= P<t>,h{Q) + E(0, 0) + «7o^(0, 0) + g 2 ^|(0, 0) + . . . 



iZq 



1 



Aq 2 + ml 



(C4) 




FIG. 11: Flow equation for the fermion self-energy. The no- 
tation is chosen as in Fig. [l] The crossed circle indicates an 
insertion of Rk in the loop integral. 



with the corresponding identification of Z, A, and m^. 
We choose Z A and A A to be nonzero at the beginning of 
the flow. In fact, they are attracted immediately towards 
a fixed point with Z k = A k = 1 in the early stages of the 
renormalization group low [ID]. Thus, choosing Z ^ 
immediately provides for the correct ultraviolet regular- 
ization of Eq. ( 12 1. Of course, we can also choose Z\ = e 



with sufficiently small e > 0. 

The fermions are treated perturbatively but in a mo- 
mentum resolved fashion. Perturbatively here means that 
we neglect the feedback of the self-energy on the other 
running couplings. Of course, they can be implemented 
iteratively, thus enhancing the quantitative precision of 
the results. We have an additional flow equation for the 
fermion self-energy shown in Fig. |11| In the second line 
of the figure we replaced the full fermion propagator by 
the microscopic one, in accordance with our perturbative 
treatment. Note that all quantities on the right hand side 
of the flow equation in Fig. [TT] are known to us and the 
fermionic self-energy can be readily integrated. 

From Eq. (Cl) we can read off the inverse propaga- 
tors for the bosons and fermions at scale fc, which are, 
respectively, given by 

h,k{Q) = + A kQ 2 /2 + ™ 2 ,k: (C5) 

P*AQ) = m + q 2 -p (C6) 

in the symmetric phase. We employ regulator functions 
R*,k{Q) = Mk 2 - q 2 m0{k 2 - q 2 /2), (C7) 



fc 2 



-(q 2 -p))9(k 2 -\q 2 -p\), 
(C8) 



which arc optimized for a derivative expansion. Note that 
the functions R k {Q) do not depend on the frequency qo, 
which facilitates a semi-analytical treatment of several 
of the loop-expressions encountered in the following. In 
addition, we introduce the abbreviations 



1 



P<p,k(q) = -jP^.kiq) = q /2 + m, pk + R<t, tk {q ), 

Pip,k(q) = q 2 - p + R^,k{q 2 )- 



(C9) 
(C10) 



We normalize the running couplings by A k and define 
S k = Z k /A k , m^ k — m^ k jA k , A^fc = X^ k /A k , po :k — 
Po,kAk, and h\ — hl/A k . Note that this introduces a 
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fc-dependence of the Yukawa coupling hi. Moreover, we 



define the anomalous dimension rjk 
substitutes the flow of A k . 



-dt log Ak, which 



Appendix D: Flow equation for the fermion 
self-energy and the contact 

The effective average action is parametrized in terms 
of the mean fields (fix and ipx, whereas the Schwinger 
functional is parametrized in terms of source fields. Both 
descriptions are equivalent and any mean field configura- 
tion, in particular with nonvanishing fermion expectation 
value, can be enforced by choosing an appropriate exter- 
nal source. Taking functional derivatives of Tk[4>, ip] with 
respect to the fields (fix and ipx serves for generating all 
one-particle irreducible correlation functions of the the- 



ory. In the same way, the flow equation ( 32 ) serves for 



generating the flow of all correlation functions. 

We consider the inverse fermion propagator for a con- 
stant (and real) bosonic background field p — <fi*<fi and 
vanishing expectation value of the fermion mean field, 



GT 1 fe (x,y,p) = 



5iPa{X) ^5^ b {Y) 



\P,0}> (Dl) 



where ipA and ifis is any of the four indices 
("01, ip2, V'l; such that we arrive at a 4 x 4-matrix. 
In a homogeneous situation, via Fourier transformation 

/ e ipx e ip ' Y G^ k (X,Y,p) = 6(P + P')G^(P,p) 

JX.Y 

(D2) 

we obtain the full inverse propagator G^ ! 1 k (P,p) in the 
presence of a background field p. 

The microscopic (or classical) inverse fermion propa- 
gator in the presence of a background field is given by 



-P t<A (-P) 
Pip,h(P) h^ifie 



(D3) 



Herein, the fully antisymmetric tensor with two indices is 
denoted by e. The physical self-energy is defined as the 
difference between the full (i.e., k = 0) inverse propagator 
and the classical one: 



= G-'{P) - GI\(P) 



(D4) 



Since for k = the bosonic field can have a nonvanishing 
expectation value po.k=o, we make this definition more 
explicit by writing 

- G^ =0 (P, p 0>fe =o) ~ G^ A (P,p , k =o). (D5) 

The most general form of the self-energy in the spin bal- 
anced case of equal chemical potentials for the hyperfine 
components can be parametrized by two complex-valued 
functions via 



E(P) 



_ /%, an (P)£ -S^(-P) 

E^(P) -(E v ,, an (P)) < 



(D6) 



Eq. (D5) is readily extended to a scale-dependent 



fermion self-energy by defining 



S fc (P) 



= G ^k( p > Po,k) - G,^ A (P, po,k)- 



(D7) 



Evaluating the self-energy for each k on the expectation 
value p ,fc 01 the bosonic field properly takes into account 
the fluctuations on different scales. As a result, the flow 
equation of the (ip*, component of E^ (i.e. the nor- 
mal contribution) is given by 



<9 t E^ fc (P) = (r]APo,k + d t po,k) ^^^ 

+ (d t \p^,k)(P,P0,k) 



(P,Po,k) 



(D8) 



with RG-time t = log(fc/A) and normalized field p — 
pAk- The derivative is performed for fixed p. On the 
right hand side, the self-energy appears as a function 
of the background field, whereas the left hand side only 
depends on k and P. 

Within our truncation we have 



/ \ 2 f R^q 2 ) (^ 2 9o.s- 2i ^o,s-2Apo(Apo+^(^)-^(g))(i(po + 9o,B)+PV'( , ? + P)) 
[d t \^, k )(P, Po ) = h,J — det 2 B (Q)det F (Q + P) 



det|j(Q)det F (Q + P) 
(iS(p +q ,F)- ^Po-P4>(q + fi)(qoF + 2[ qo,FP4,(q) + h 2 p -p'l(q)) 

RA<f)- : — - (D9) 

det B (Q + P)dct F (Q) 



with g ,s = 27rnT and go,F = 2n(n + 1/2)T being 
bosonic or fermionic Matsubara frequencies, respectively, 



and po = 2ir(m + 1/2)T. We introduced 

dets(Q) = S 2 ql J3 +P4,(q)(p4,(q) + 2Ap ), 
det F (Q) = ql.p +P%{q) + h^po 



(D10) 
(Dll) 
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and adopted the notation of Eqs. ( C9 ) and ( CIO ) for p^^ 
and p^.k. The Matsubara summations can be performed 
analytically and we arrive at an explicit expression for 
the third term in the flow equation ( |D8[ ). 

The flow equation for E^,fc(-P) is valid for arbitrary 
values of P. For large P, the equation simplifies consider- 
ably and, eventually, allows to derive the renormalization 
group flow of the contact Ck- We restrict the following 
discussion to the first integral in Eq. (D9|, which is re- 



sponsible for the high momentum behavior. We have 



i(po + go)+MQ + P) 
detjr (Q + P) 

for large p 2 and find 



1 



'ip,k 



with 



{dt\pc k )(p) = h\ 



A 



-ipo + p 2 - (i 



(d t \pc k )(p) 



(D12) 



(D13) 



s2 iIm - 2[ Sq 0iB - 2Xp (Xp + p^(q)) -p$,(q) 



det|(Q) 



i> 2 i 'fry. 



2S 2 u 2 



4>,k 



2X 2 p 2 



detfl(Q) 



de4(Q) 



For the definition of w^fc see Eq. (44) 

ntact ac 

-Ck{po,k 



We define the flowing contact according to 
1 

T 



Ck 



(D14) 



(D15) 



where Ck{p) '■= Ck(p) is expressed in terms of the normal- 
ized field p — pAk- Since 



c(k,p) = c(k,p(p, k)) 

we have 

dc dc dc 

dk dk dp 

dc 1 dc 

dp A k dp ' 



dp 
dk 



dc 

dk + k 



1 _dc 

VAPl 



dp 



(D16) 

(D17) 
(D18) 



From these two relations we deduce the flow equation 
for Ck (p) in the presence of the fc-dependent background 
field p = Akp to be given by 



uh\„n n,>) = «H/v. >i/> = />t/.)) + >iAP^{p). (mi)i 



Thus we arrive at 
d t Ck ' 



! y(dt\pCk)(po,k) + (vaP + d t pQ,k)-g^(po,k) 



(D20) 



We show that <4(po,fc) = h 2 hk to leading order in E^. 
Indeed, the fermion propagator in the presence of a real 
background field within our truncation is given by 



G^{Q,p) 



hcf>e P,p,a{-Q) 

detjr(Q) y-P^,AQ) -He 



with 



det F (Q) = PtpAQ)PipA-Q) + h2 P- 



(D21) 
(D22) 



The following relations are valid up to (perturbative) 
corrections due to the fermion self energy, i.e. to order 
O(S^). We find 



dG, 



(Q) 



h/2(j) 







dp det F (Q)" v ^ / ' detjr(Q) V° ~ £ 

and, accordingly, for the inverse propagator 

d ^(Q)-~G^(Q) d ^(Q)G^(Q) 

_ h l r -i (n] W A o 

- d^Q) G * (Q) ~ 



(D23) 



det F (Q) V° 1 



With 



G ^ (Q) - \P*pAQ) 



h(j>e -PfA~Q) 

h<pe 



(D24) 



(D25) 



we then find for the , ipi )-component at large external 
momenta \Q\ ^> h 2 p the relation 



dE^ 

dp 



(p) 



P+AQ) p i>A-Q) 



p i>AQ) 



h 2 



p^A-Q) 

(D26) 



This completes the proof of 
dc 



dp 



(p 0tk ) = h^ + O(^). 



(D27) 



The flow equation for the contact becomes 
hi 



dtC k ■ 

4 



-(ijAPo.k + d t pQ,k) 



Rd 



i 



2X 2 p 2 



dets(Q) 



dct^(Q) 



(D28) 

The result of the Matsubara summation and the q— 



integration is given in Eq. ( 43 1 



To clarify the relation between dtCk and dtn^fa we 
consider the scale-dependent density n& defined by 



»k = 



dUk 

dp 



(Po,k)- 



(D29) 
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For k — we arrive at the physical density 
dP(p,T) dU(p,T,p ) 



dp 



dp 



(D30) 



with pressure P{p,T). The corresponding flow equation 
is given by 

d t n k = -a k (riAPa,k + d t p Q ^) - (dndt\pUk\{pa,k)- 

(D31) 



with 



a k = o a \P0,k)- 



dfidp 



(D32) 



The flow of the effective potential receives contributions 
from bosonic and fermionic fluctuations. Defining 

= -( d M U i B) ) M, (D33) 



°t n ^,k 



■t\f> u k 



(p 



•0,k) 



(D34) 



as the bosonic or fermionic contribution, respectively, we 
can employ 



(dt\pUi B) ){p^ 



h 

Rd 



k- 2 + U' k {p)+ P U' k \p) 



q A [k2 + U>(p)+2pUZ(p)][k 2 + U>(p)}+S*q 2 

(D35) 



to find 



d t n <b,k = ~ a k 



f &±( I 

Jo A Vdets 



2S^l k +2\ 2 p% 



(Q) 



det|(Q) 



(D36) 



This is precisely the third term in the flow equation (D14) 



for Cfc. We can summarize these findings in the generally 
valid relation 



d t c k = *-(-~. la-.-^xaJ") 



4a fe V 



-ak(VAPo,k + d t p ,k) + dtn^J, 



(d t n k - d t n:Yl) 



4>,k J i 



(D37) 



where n 



(U) 

<t>/ip,k 



are contributions which arise from effects 



of bosonic/fermionic fluctuations on the effective po- 
tential and hence the density. In contrast, the term 
—o: k (riAPo,k+dtPo,k) accounts for the nontrivial scaling of 
the renormalizcd propagator and the contribution from 
condensed bosons. 



the flow of Cfc is zero, because /?o,fc = and there are no 
bosonic fluctuations on high energy scales. Hence 



d k C k = forfc>/i 1 / 2 ,T 1 / 2 



(D38) 



However, the flow of the prefactor h^/at = h\/a k is gov- 
erned by the scale set by the scattering length a -1 . On 
the far BEC side, this quantity is large and these renor- 
malization effects set in far above the many-body scales. 
At such high scales the vacuum relation a k — —2Zfak, 
with wave function renormalization Z$ k of the bosons, 
is valid. The relation stems from the appearance of the 
combination Z r j >tk {d T — 2p) in the propagator due to semi- 
local U(l)-invariance and the symmetry preserving na- 
ture of the flow equation. One can show that the vacuum 
flow of Zfck is solved by Z$,k=o = h^a/32ir on the BEC 
side 38 . i.e. 



4a, 



-> — forfc» M 1 / 2 ,T 1 / 2 . 



(D39) 



This prefactor effectively enters the renormalization 
group equation of C k which takes place on the many- 
body scales and thus we can write 



dkC k 



-in 



(D40) 



Due to the fact that there are no fermion fluctuations 
contributing to the density on the BEC side, we have 
n^_k = n k and thus arrive at the well-known relation 
Cbec = 47rn/a. The mean- field result Cbec — 47m/ a 
receives corrections from bosonic fluctuations, which are 
incorporated in the renormalization group flow. 



Appendix E: Perturbative contact from the equation 
of state 

We consider the perturbative BEC regime of the 
crossover, i.e. the region 



(k F a) 1 > 1. 



(El) 



Moreover, we restrict the considerations to the case of 
zero temperature. We closely follow Tan's presentation 
in Ref. [8]. The energy density of the system (in the 
canonical variables) is given by the Lee-Huang- Yang ex- 
pression 



E 
V 



n d e B 



-nAl 



128 



n d a\ 



(E2) 



where n<j is the number of dimer atoms, which is related 
to the density according to na = n/2 on the far BEC 
side , and Eb is the binding energy of a dimer. From Eq. 
(El) we find the gas parameter n^a^ to be small. The 



In the BEC limit, Eq. (|D37j) simplifies considerably. coupling constant for the dimers is given by 



Indeed, following the flow of C k from k = A to k = 
we see that in the early stages of the flow, where k is 
much larger than the many body-scales set by p and T, 



9d 



(E3) 
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where we used that Md = 2M = 1 and = na. The di- 
mensionless constant k relates the scattering length of 
dimers to the fermionic scattering length. The exact 
value of k is known to be 0.6 [58) . However, within our 
truncation we have k = 0.72 from a solution of the vac- 
uum problem [40]. We employ the latter value for con- 
sistency and emphasize that k is not a free parameter in 
our model. The energy density in terms of n and a is 
found to be 



E 
V 



SB 



1 



128 



:V nn 3 a? 



(E4) 



We compute the chemical potential according to fi(n) — 
d(E/V)/dn and subtract half the binding energy to ob- 
tain the (positive) many-body chemical potential 



Aimb(n) =M- y = 71 nan 



fl 



32 
3V2tt 



V nn 3 a 3 



(E5) 



Inverting this relation to the same order of approximation 
we find 



mb ) 



irna V 



32k 



(E6) 



This is the equation of state in the grand canonical vari- 
ables at T = to leading order in the gas parameter 

V^mbO. 

By virtue of the adiabatic sweep theorem ^ we can 
derive an expression for the contact in the BEC regime. 
For this purpose we employ £b = — l/Ma 2 = — 2/a 2 . 
From Eq. ( E4 1 we then find 

Cbec 



inn 



+ K7r 2 a 2 n 2 ( 1 + 



G4 



Vr. 



•■3/^3 



3V27T 

at zero temperature. By dividing by k F we obtain 
Cbec 



(E7) 



3tt 



(kpa)~ 



i 



64k 3 / 2 



37rv67r 



(kpa) 



3/2 



Inserting the equation of state ( E6 ) we arrive at 
Cbec 



u 2 

AVb 



-{yfjhnba) 2 f 

K V 



1 - 



32k 
3tiV2 



'MmbO 



(E8) 



(E9) 



We now consider the asymptotic behavior on the BCS 
side. The equation of state at zero temperature is given 
by 



E 3 2 / 
= -kwn[ 1 



10 



V 5 r V 9tt 
We deduce for the chemical potential that 

4 



d E 
dnV 



K 1 



3tt 



-kva 



(E10) 



(Ell) 



The contact is found from 



d(-l/a) 



(E12) 



Thus we have 



Cbcs 
hp 



9tt 2 



Inserting the equation of state k-p (/i) 
at 



Cbcs 



9tt 



(fc F a) 2 . (E13) 
fjl we then arrive 

-(^a) 2 . (E14) 
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